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I. INTRODUCTION 

Classical chaos as an unpredictable motion in systems with causal dynamics is 
well understood phenomenon |]T]-§[. Many approaches have been suggested during 
the past decade to define chaos in quantum systems. Nevertheless, this problem is 
far from clear [[||J. We will use in this work the terminology of Refs. [|5]J|, where 
quantum chaos has been defined as the study of semiclassical behaviour of parameters 
of quantum system whose classical analog exhibits chaos. 

Very interesting models in the study of quantum chaos properties are models of 
quantum optics, among which are most important ones describing interaction between 
atoms and light [[/J. The Jaynes - Cummings model (JCM) || of a single two - 
level atom interacting with a single - mode cavity photon field attracted widespread 
attention of many investigators p|-|TTf. It has been observed theoretically fT2| , |T3| that 



the model possesses a considerable and unexpected pure quantum features (periodic 
collapses and revivals of the mean atomic exitation energy and the mean photon 
numbers). 

Belobrov, Zaslavsky and Tartakovsky |14J first studied the dynamics of the many 
- atom JCM at the semiclassical approximation. Later several variants of this model 



have been investigated by different researchers Hl5|,[7|, p!6| , p!7| . It was found that the 
chaos came into view in this system when the rotating wave approximation (RWA) is 
not used, and this implies the very large values of transition dipole moments and atom 



densities. Gorokhov and Ruchkov Jl8| was observed that semiclassical chaos may be 
efficiently suppressed when possible parametric instabilities of nonlinear media packed 
in the cavity are taken into account. 

The possibility of chaotic behaviour for system which consists of the single two- 
level atom coupling with the one - mode of quantized electromagnetical field in ideal 
resonator, is also discussed at the last time [|llj . Contrary to the many - atom systems 
in this situation the ratio of atom - field coupling constant to transition frequency is 
usually very small and chaotic regime cannot appears. 



Up-to-date technological achievments [§3J lead to practical interest for systems 
with coupling constant of order of frequency u, when RWA can not be applied. 



The experiments with the one-atom maser [20, 21 1 and recent realization of an 



optical microlaser P2| , |23| gives hope to achieve more favourable values of this ratio 
and to observe chaos for this fundamental quantum - electrodynamical system. 

Besides, it is interesting to consider different quantum optical models with the aim 
of checking their possible quantum chaotic properties. The two - level atom with m— 
quantum transitions should be studied first. Another case is a quantum model of the 



three - wave parametric interaction [24]. As we show here, it leads to hyperbolic JCM 
- a system with model hamiltonian in terms of SU(1, 1) group generators instead of 
the Pauli matrices describing two - level atom. 

The paper is organized as follows. Section 2 presents a brief description of model 
hamiltonians and basic assumptions. Section 3 is devoted to the introduction of 
relevant group theoretical coherent states (CS) and their applications to the derivation 
of motion equations using CS path integral representation for a transition amplitude. 
The most important feature of coherent states is thatthey are pure quantum states, 



which are the most close for classical ones [25,26] . The coherent-state point's motion 



obeys classical-like dynamical equations, and it is reasonable to use CS basis for 
invoking chaotic mystery into quantum mechanics. Section 4 deals with the results 
of computer simulation of systems discussed. In Conclusions we have considered 
possible generalizations of model hamiltonians and developing approaches for purpose 
of taking into account quantum corrections and their influence on chaos. 

II. MODEL HAMILTONIANS 

Let us consider the two - level atom model and some of its generalizations. It is 
well known, that the quantum hamiltonian for two-level atom, coupling one mode of 
quantized electromagnetical field in an ideal cavity, can be represented as the following 
0: 



H = 


= H F + H A + ( ~H. int + ~Hi nt ) , 


H F = 


= ,(b+b + i), 


H A - 


= u3 , 


tt(1) 


= ghJ + + gb + J + + h.c. , 


tt(2) _ 
"■int ~ 


= A;(b+ + b) 2 , 



(2.1) 
(2.2) 
(2.3) 
(2.4) 
(2.5) 

here v is (angular) frequency of the field, uj is atomic frequency and difference 
between its energy levels, as we choose % — 1, g and k are coupling constants, 
b + and b are the boson creation and annihilation operators, Jo, J± = Ji ± J2 
are the generators of SU{2)— group, which equal one half of corresponding Pauli 
matrices (only two levels are taken into account and operators J± increase (decrease) 
the energy of the atom by u ) . 

In agreement with the model, v w u and the resonant terms in the Hamiltonian 

Hi„£ gi ye principal contribution, when amplitude of transitions is calculated in the 

(2) 
case of low values of constant g. Usually, the term H^ is eliminated, because k -C g. 

Also, g b + J + and hermitian conjugate term discard in rotating wave approximation 

(RWA). Finally, Jaynes - Cummings Model (JCM) is come to. 

Notice that the antiresonant terms g (b + J + + bJ_) have appeared in similar 

hamiltonian for electron - phonon interaction in solids. The Jaynes - Cummings 



Model without RWA is called in the literature [|7| the "dressed" JCM. 

As an evident generalization of " dressed" JCM, we may consider first the hamil- 
tonian of a (2 j + 1)— level atom (spin j = 1/2, 1, . . .) in the form of Eq.(2.1) with 
the (2 j + 1)— dimensional generators of SU(2) group. Secondly, in many cases m— 
quantum transitions between atomic energy levels are of present interest. It leads to 
nonlinear in bosonic operators generalizations of interaction hamiltonian 

H« = g (b) m J + + g (b+) m J + + h.c. (2.6) 

It is important to notice for the future considerations that all these hamiltonians 
can be expressed as an operator valued functions of the direct product of Lie groups 



W(l) ® SU(2) generators (the first group corresponds to field (photon) subsystem 
and the second - to atomic subsystem.) 

The atomic subsystem's dynamical symmetry group 577(2) has the non-compact 
group 577(1,1) as an analitic continuation p6| . Unlike compact 577(2) with unit 



sphere as quotient space, quotient space of SU(1, 1) is a hyperboloid. The unit cir- 
cle's interior is suitable substitution for this space at complex plane, its boundary 
meets Lobachevsky plane infinity Hamiltonian with generators SU(1, 1) instead of 
generators SU(2) is appearing to describe three-wave parametric interaction in non- 
linear media. The main difference between SU(2) and SU (1,1) ("hyperbolic" JC) 
models consist here in the role of nonresonant terms. Their presence in hamilto- 
nian actuates a parametric build-up of oscillations and leads to energy and number 
of quanta increase in all subsystems. Therefore, the process requires external time 
dependent sources. 

It is cause to write down Hamiltonian of the model as the following 

H = z/(b+b + -) + 2 uK Q + g K+b + gK h + + AK+b + AK b+, (2.7) 

K and K-t are generators of SU(1, 1)— group, ( [K , K±] = ±K±, 
[K + ,K_] = — 2K ); A is (the strength of the external field - system interaction) 
determined by properties of external classical electromagnetic field and v ~ 2u is 
treated as a condition of resonance. (Generally, 2u = uj\ +002, where uj\ t 2 are the fre- 
quencies of electromagnetic waves parametrically coupled with wave of the frequency 
v, and constant g is proportional to some component of tensor of the third order 
nonlinear susceptibility.) 

Notice, here it is not necessary for A be equal of g, and in general, A is not a 
constant. 

Considering that all introduced model hamiltonians possess dynamical symmetry 
group described above, coherent state technique would be applied. The technique for 
linear of dynamical group generators hamiltonian permits us to get explicit equations 
for coherent state parameters' motion, unlike in our case approximate equations are 



expected only. In our approach coherent state means the direct product of group - the- 
oretical coherent states of W(l), and "atomic" dynamical group (577(2) or SU(1, 1)) : 

\CS>= |a>®|C>, (2.8) 



here \a > is usual Glauber CS and \( > is "atomic" CS [26]. 



III. COHERENT STATES, PATH INTEGRALS AND EQUATIONS OF 

MOTION 

The purpose of this section is to get a semiclassical equations of motion in coherent 
state representation for JCM and its generalizations. Most natural way to do that 
consists in using a path integral formulation of the problem [p^ ,^]. 

Consider a quantum system with the hamiltonian being an operator - valued 
function of the generators of the dynamical group W(l) ® 577(2) (577(1, 1)) and use 
Schweber's Method. Here it is described briefly. Introduce a matrix element of the 
evolution operator of a system U(t, t ) = exp(— z(£ — £ )H), between group - theoretical 
coherent states 

< CS\exp{-i{t - * )H|C5' >= 

= [■■■[< CS N \ expH (* - *o) H/iV)|C5 iV _i > dfi(CS N ^)x 

(N-l) 

< CS N -i\exp(-i(t - t )n/N)\CS N ^ 2 > d^{CS N ^ 2 ) x ... 
... x dfi(CSi) < CSt\ exp(-i(i - t )H/A^)|C5 >, (3.1) 

At = ^, iV»l, 

d/j,(CSi) is the invariant measure in the resolution of unity. Here group property of 
evolution operator and completeness condition are involved. Each factor in Eq.( |3.1| ) 
is transforming: 

< CS j+ i \e- t{t - ta) ™ /N \CS 3 >~< CS j+1 11 - t{t - to)H/N\CSj >~ 



~< CS j+1 \CSj > exp(-i(t-to)h(CS j+1 \CSj)) 

h(KC j+1 \KCj) 



< CSj+i\H.\CSj > , , ... ^ 

<cs J+1 \c Sj > = h (°^l°-0 



is a covariant symbol of Hamiltonian p9 



Then we go over to limit by N — > oo and have as result: 



< a,(\exp{-i{t-t )H)\a/,( f >= / DM^e 



S - is treated as action function, and Euler - Lagrange equations take the following 
form after some algebra for system with symmetry group W(l) <S> 577(2) : 



a + i va + i 2jgm(a*) m ~ l T ^ + i 2j\m(a*) m ~ 1 T ^ = 



a* - i ua* - i 2jgm{a) m - 1 T £^ - i 2jAm(a) m " 1 I ^ = 



i+ee 



(3.2) 



C + iuC + i g{a) m + i A(a*) m - i g(a*) m ( 2 - i \(a) m ( 2 = 



C* -iuC - i g(a*) m - i \(a) m + i g(a) m C 2 + % A(a*) m C* 2 = 0, 
here (2 j + 1) is the number of atomic levels. 

Note, so symbols * as ~ designate complex conjugation, the first for coherent-state 
parametres and the second for coupling costants. In the following calculations we 
are considering two particular cases: g ^ 0, A = - (JCM), and g = A - ("dressed" 
JCM). 

Similar equations are led to system with dynamical group W(l) ® SU(1, 1) : 

a + iva + i 2kgm(a*) m ~ l j^* + i 2k\m(a* 



\m— 1 C* _ n 
I 1-CC* U 



a* — iva — i2kgm(a) m 1 jzf^ — i2k\m(a) m 1 ^^ = 



C + icuC + i g{a) m + i \(a*) m + ig(a*) m ( 2 + i \(a) m ( 2 = 



C* - i wC* - ig(u*) m - % \{a) m - i g{a) m C 2 - i X(a*) m C 2 = 



(3.3) 



Instead of j here arises the parameter k related with eigenvalue of the 
577(1, 1)— group Casimir operator K 2 = K 2 , — 1/2 ( K + K_ + K K + ). For two-mode 
oscillator k = 2 = 1/2, 1, 3/2, . . . (An is the difference between the number of 
quanta in modes); in one- mode oscillator with two-quantum transitions two possible 
values exist: k = 1/4 for even levels, and 3/4 for odd levels. 

These equations are non-linear ones and have explicit analytical solutions in some 
trivial cases only. Therefore, they should be solved by computer calculations. 

Using Dormand - Prince Method of fifth order computer program permits us to 
get graphics of real and imaginary components of coherent states parametres' time 
evolution, complex - plane evolution of CS - parametres both for field (in a plane 
and for "atom" (in the complex plane realizations of quotient spases of dynamical 
groups SU{2) and SU(1, 1). Evolution of atomic (S77(2)) coherent state on the 
Bloch-sphere is also realized when the complex parameter ( is replaced by two angle 
parameters with the substitution: ( = e %< ^ cot(|) that gives a stereographic map 
of a complex plane to unit sphere |26| with the South Pole corresponding to the 



origin of the plane. Moreover, phase portraits of real and imaginary coherent state 
parameter components and the transition probability time-dependence are received. 
For definiteness we restrict in the paper the consideration of the case to find atom in 
"down" state: 

P{t) = | < down\({t) > | 2 . (3.4) 

IV. RESULTS OF COMPUTER SIMULATIONS 

In present section the results of computer calculations are described. 
It is well known, Jaynes - Cummings Model has regular dynamics, for example, 
up-down transition probability's behaviour obeys formula |T0| 



q 2 (n+l) . 2 



(is - lu) 2 + q 2 (n + 1)- _ y 



-fi (t-to) 



(4.1) 



if number of photons in the initial state is fixed. Analogous but more unwieldy 
formulae exist for some other type initial states, including coherent state. 

However, calculation with non-resonance terms leads to considerable complication 
of time-dependence behaviour so transitions probability as coherent-state parameters. 
This complication is most conspicuous in case of a small average number of photons 
(< n >~ 1) in the initial state of field mode and large coupling constants. The 
following regularity might be noted: while g < 0.05 -f- 0.1 (magnitude of coupling 
constant depends on initial values of atomic and field coherent - state parameters) 
transitions probability time behaviour remains "sine-similar" curve, then with in- 
creasing g "sine-similar" envelope of curve exists only and multifrequency regime is 
observed. Further, if coupling constant is continuing its growth, there is no harmon- 
ical envelope of the curve, but such value of g exists when number of probability 
oscillations frequency cuts down (to two frequencies, in our case) (as Fig.l, which 
has been calculated with g = X = 0.5, v — u — 1, , and system has been prepared 
in the initial state: |a(0) > <8>|C(0) > , C(0) = 0.3, a(0) = 1.). After this, system 
(in sense of probability behaviour) returns to multifrequency regime (Figs. 2, 3) and 
draws near the chaotic regime (Fig. 4). The effect takes place with < n > essentially 
more then 1 (about 25, in case is considered), what have been described. The effect 
is observed on lesser value of the coefntient g then it is necessary for small average 
number of photons in the initial field mode. It is important then < n >^> 1 corre- 
sponds to classical field and the problem reduces itself to problem of "spin dynamics 
in the (magnetical) field" , with no chaotic solutions. The field coherent state point is 
interested to move over constant radius circle with very small distorsions. Also, case 
of two - quantum (2u — u — 1) transitions "dressed" JCM have been investigated. 
Some typical results are showed in Figs. 5 - 9 (here g = X = 0.4, £(0) = 0.3). The 
first picture (cc(0) = 1) of them presents chaotic behaviour of system unlike second 
one with larger initial field coherent state value («(0) = 5) is multifrequence regime. 
These conclusions are founded by observing of calculating phase portraits (see Figs. 7, 



8). The last picture (Fig. 9) in series demonstrates the beginning evolution coherent 
state point over Bloch - sphere. 

In conclusion, systems behaviour in both described cases are similar, but in two- 
quantum "dressed" JCM a chaotic regime arises with in order smaller of coupling 
constant magnitude. Although, we guess that needed value achievement in real ex- 
periment is very difficult problem. 

The possibility of "dressed" JCM chaotic behaviour is provided by antiresonant 
terms, which might describe as single virtual two-photon process. Similar terms may 
be included in hamiltonian of three-wave parametric interaction, but their source 
must be found out of system, because energy conservation law (for closed system) 
is violated by the terms. The medium with non-linearity of forth order is occurred 
the cause of parametric build-up process, for example. The short-time exposure is 
simulated by step-function A(t) : this function identically equals zero everywhere 
except short (as compared with evolution time) interval. It is occured, this way does 
not lead to chaos. The changes in phase portraits are observed, and there is a marked 
difference between some of portrait, but all of them stay as regular structures (Figs. 
10 - 13). The atomic coherent state point approaches to bound of unit disc for the 
finite time, and the time decreases with increasing coupling constant or interaction 
time interval. 

The same effects are observed for time - dependence of coherent state variables. 
For example, the Figs. 14, 15 show dependence Re(, when activating the stepwise 
parametric interaction with following parameters of system: v — 2oo — 1, g — 0.2, for 
initial state determined by £(0) = 0.5, a(0) = 5; the interaction was switched after 
ooto — 10 free evolution of system for interval S(uto) = 0.2, A = 0.2. 

The destruction of system is a possible result of exposure parameters increasing. 
Physical treatment of this is exponential growth of quanta average number in so field 
as "atomic" modes. As result, hyperbolic JCM is simpler in "chaotic sense". Per- 
haps, the behaviour of system nearly parametric generation domain requires detailed 
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studying by more accurated computation procedure unlike used. 

V. CONCLUSIONS 

In the paper we restrict our discussion only particular cases of the models (two 
- level atom (j = 1/2) or degenerate "atomic" system with dynamical symmetry 
group SU(1, 1) (k = 1/4 or k = 3/4)). The computer program allows to do similar 
calculations for wide class of possible evident model enlargement. 

Besides, the developed approach gives a natural way to its generalization for 
chaotic quantum corrections' calculations, and presented results might be consider 
as the first step of our investigation of chaos in quantum optical models. 

Method of Heisenberg operator averages, based on similar ideas (but non tech- 
nique), gives the opportunity to involve wider class of problems because it don't needs 
of assumption about conservation subsystems' states as coherent during the evolu- 
tion. Nevertheless, account of quantum corrections is believed to be more complicated 
technically in this approach. 

We compared the results in both approaches for "dressed" Jaynes - Cummings 
Model by calculating inverted population difference. All of them are in good agree- 
ment. Comparison of these methods and computer simulations for some extension of 
quantum optical model hamiltonians will be published in the next paper. 
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Figure Captions 

Fig.l. Time dependence of the down - state population in two - level atom 
with g = A = 0.5, v = u = 1, C(0) = 0.3, a(0) = 1. 

Fig. 2. Time dependence of the down - state population in two - level atom 
with g = A = 0.6, v = u = 1, C(0) = 0.3, a(0) = 1. 

Fig. 3. Time dependence of the down - state population in two - level atom with 
g = A = 0.7, v = w = l, C(0) = 0.3, a(0) = 1. 

Fig. 4. Time dependence of the down - state population in two - level atom with 
g = A = 0.8, v = u = 1, C(0) = 0.3, a(0) = 1. 

Fig. 5. Time dependence of the down - state population with the two-photon 
transition in two - level atom and g = A = 0.4, 2 v — to — 1, ((0) = 0.3, a(0) = 1. 

Fig. 6. Time dependence of the down - state population with the two-photon 
transition in two - level atom and g = A = 0.4, 2 v — u — 1, £(0) = 0.3, a(0) = 5. 

Fig. 7. Phase portrait in the plane (Re (; Re( ) (g — X — 0.4, 2 v = uj = 1, a(0) = 
1,C(0) = 0.3.) 

Fig. 8. Phase portrait in plane (Re(; Re() (g = A = 0.4, 2 z/ = u; = l,a(0) = 
5,C(0) = 0.3). 

Fig. 9. Bloch - sphere atomic coherent state evolution (g = A = 0.4, 2u = u = 
l,a(0) = l,C(0) = 0.3). 

Fig. 10. SU(1, 1)— coherent state dynamics on unit disc \(\ < 1, when activating 
the stepwise parametric interaction (k = 1/4, v = 2u = l,g = 0.2, a(0) = 1,C(0) = 
0.5; <$(urf ) = 2, cut = 10, A = 0.5). 

Fig. 11. Phase portrait in plane (Re(; Re() when activating the stepwise para- 
metric interaction (k = 1/4, v — 2uj — l,g — 0.2, a(0) = 1,C(0) = 0.5; S(ut ) = 
2, ujt = 10, A = 0.5). 

Fig. 12. SU(1, 1)— coherent state dynamics on unit disc \(\ < 1, when activating 
the stepwise parametric interaction (k = 1/4, v — 2uj — l,g — 0.2, a(0) = 1,C(0) = 
0.5; S(ut ) = 3, ut = 10, A = 0.5). 
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Fig. 13. Phase portrait in plane (Re(; Re() when activating the stepwise para- 
metric interaction (k = 1/4, v — 2uj — l,g — 0.2, a(0) = 1,C(0) = 0.5; S(uto) = 
3, ut = 10, A = 0.5). 

Fig. 14. Time dependence Re(, when activating the stepwise parametric inter- 
action (k = 1/4, v = 2 uo = 1, g = 0.2, a(0) = 5, C(0) = 0.5; S(ut ) = 0.2, ujt = 
10, A = 0.2). 

Fig. 15. Time dependence Re(, when activating the stepwise parametric interac- 
tion (k = 1/4, v = 2u = l,g = 0.2, a(0) = 5,C(0) = 0.5; S(ut ) = 0.2, uot Q = 10, A = 
0.2.) (The initial part of Fig.14., out < 40.0). 
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